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Abstract. Models of steady-state plane-parallel shock 
waves propagating through the unperturbed hydrogen 
gas of temperature Ti — 6000K and density pi — 
10"^" gm cm~^ are computed for upstream velocities 
15 km s~^ < C^i < 70 km s~^. The properties of the ambi- 
ent gas are typical for atmospheres of pulsating stars. The 
shock wave structure is considered in terms of the self- 
consistent solution of the radiation transfer, fluid dynam- 
ics and rate equations for 2 < i < 4 atomic bound levels 
with a continuum. The radiative flux emergent from 
the shock wave was found to be independent of the lower 
limit vl of the frequency range provided that vl < 1^2, 
where 1/2 is the Balmer continuum head frequency. At the 
same time the decrease of vl is accompanied by decrease 
of the Lyman continuum flux and leads to smaller heating 
and weaker ionization of the hydrogen gas in the radiative 
precursor. For all models the size of the radiative pre- 
cursor is of ^ 10** cm and corresponds to several mean 
free paths of photons at the frequency of the Lyman con- 
tinuum edge vi. The compression ratio at the discontin- 
uous jump gradually increases with increasing upstream 
velocity Ui, reaches the maximum of p'^ / p~ = 3.62 at 
[/i « 55 km s~^ and slowly decreases for larger Ui due to 
the strong rize of the preshock gas temperature. The radia- 
tive flux from the shock wave was determined as a function 
of the upstream velocity Ui and its ratio to the total en- 
ergy flux in the shock wave 62 was found to range within 
0.18 < Fr/Cz < 0.92 for 15 km s'^ < Ui < 65 km s'^. 
Thus, at upstream velocities Ui > 60 km s^^ the shock 
wave losses more than 90% of its total energy due to ra- 
diation. For all shock wave models the role of coUisional 
processes in both bound-bound and bound-free atomic 
transitions was found to be negligible in comparison with 
corresponding radiative processes. 

Key words: Shock waves - Hydrodynamics - Radiative 
Transfer - Stellar Atmospheres 



1. Introduction 

Attempts to obtain the self-consistent solution of the 
equations of fluid dynamics and radiative transfer for 
shock waves propagating through atmospheric gases were 
undertaken during several last decades. The need to em- 
ploy the self-consistent solution is of tremendous impor- 
tance in the case of the partially ionized preshock hydro- 
gen gas undergoing the substantial radiative heating and 
photoionization due to the strong absorption of the Ly- 
man continuum radiation outgoing from the postshock re- 
gion. In order to correctly treat the coupling between the 
preshock and postshock regions the two different meth- 
ods known as the "mean-photon approximation" and the 
"asymptotic layering method" were previously proposed 
(see, for recent review, Gillet 1999). Unfortunately, in both 
these methods the radiative transfer is treated as an ini- 
tial value problem which is solved with use of the shooting 
method. The principal disadvantage of such an approach 
is that the true solution of the transfer equation is un- 
stable for optical depths r > 1 since it is affected by the 
exponentially growing errors. This obstacle, however, can 
be overcome if the radiative transfer is solved as a two- 
point boundary value problem. In order to reach the self- 
consistency between the gas flow and the radiation field 
Fadeyev & Gillet (1998, hereinafter referred to as Paper I) 
employed the procedure of global iterations. Each cycle 
of the global iteration procedure involves the solution of 
the initial-value problem for the fluid dynamics and rate 
equations which is followed by the solution of the two- 
point boundary value problem for the radiation transfer 
equation. As was shown in Paper I, the global iterations 
converge and provide with the stable self-consistent solu- 
tion. 

The goal of this study has been to further develop a 
method of global iterations for calculating the structure of 
radiative shock waves. In Paper I the study was restricted 
to the model of the hydrogen atom with two bound levels 
and a continuum. The rate equations were solved for the 
first bound level and free electrons, whereas the second 
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atomic level was treated in LTE. In our previous study 
we did not take into account also the rise of the elec- 
tron temperature due to the adiabatic compression at the 
discontinuous jump. Moreover, the convergence of global 
iterations has been reached for shock models with insuf- 
ficiently thick postshock region. As a result, the models 
presented in Paper I did not cover the full recombination 
of the hydrogen and estimates of the total radiative flux 
emergent from the shock wave remained uncertain. 

In the present work the system of the ordinary differ- 
ential equations is extended by the terms taking into ac- 
count the contribution of the radiation field in the energy 
and momentum equations. The rate equations are solved 
for all bound levels of the adopted hydrogen atom model. 
The substantial improvements of the computer program 
allowed us to obtain the good convergence of global iter- 
ations for shock wave models involving almost the entire 
recombination zone and, thereby, to obtain the reliable es- 
timates of the radiative flux emergent from shock waves 
with upstream Mach numbers as large as Mi « 7. As 
in Paper I, the present study is confined to shock waves 
propagating through the partially ionized hydrogen gas. 

The present paper is organized as follows. In Sect. ||we 
describe details of our shock wave model, give the fluid dy- 
namics and rate equations written in the form of the ordi- 
nary differential equations and describe the solution of the 
transfer equation for the both continuum and spectral line 
radiation. In comparison with Paper I the equations are 
completely rewriten with the use of more convenient nota- 
tions and should replace those of Paper I where numerous 
typographical errors resulted from the missing of the TJt;X 
macro definition file. In Sect. |^ we discuss the properties 
of the solution, demonstrate its convergence and degree 
of approximation. General properties of the shock wave 
models are presented in Sect. ^. More detailed description 
of the models is given in Sects. ^ - where we discuss the 
properties of the radiation field and the structure of the 
preshock and postshock regions, respectively. 

2. Basic equations 

We consider steady-state, plane-parallel shock waves 
propagating through the homogeneous medium composed 
solely of an atomic hydrogen gas. The solution is deter- 
mined as a function of three parameters: the unperturbed 
gas temperature Ti, the unperturbed gas density pi and 
the velocity Ui at which the shock wave moves with re- 
spect to the ambient gas. The model is represented by 
a flat flnite slab comoving with the shock wave and the 
origin is set at the infinitesimaly thin discontinuous jump 
(adiabatic shock front) where thermodynamic quantities 
undergo an abrupt change. 

As was shown in Paper I, the stable self-consistent 
solution for the radiative shock wave structure can be ob- 
tained with the method of global iterations. In this method 
the rate equations for atomic level populations and equa- 



tions of fiuid dynamics are treated as an initial-value 
problem for the system of ordinary differential equations, 
whereas the radiation transfer is solved in the framework 
of the two-point boundary value problem. The interplay 
between the gas flow and the radiation field is taken into 
account by the iteration procedure. Such an approach al- 
lows us to substantially diminish the role of the stiffness of 
the rate equations as well as to avoid large errors arising 
due to the strong frequency dependence of the radiation 
field. 

In order to apply the Feautrier technique (Feautrier 
1964) for solution of the transfer equation the slab is di- 
vided into iV - 1 cells (see, e.g., Mihalas 1978). The cell 
sizes increase in both directions from the discontinuous 
jump locating at the 3-th cell interface with space coor- 
dinate Xg = 0. In cell interfaces with space coordinates 
Xj, where 1 < J < -/V, we define the monochromatic opti- 
cal depth T,jj and the monochromatic radiative flux F,jj at 
frequency v. Other variables are deflned at the cell centers 
with space coordinates 

^,-1/2 = ^(^,-1+^,), (j = 2,... ,iV). (1) 

The outer boundaries of the preshock and postshock re- 
gions are Xi and X^, respectively. For these boundary 
points we write 

^1/2 — ^1, — (2) 

since all variables must be also defined at both boundaries. 

Integration of the ordinary differential equations be- 
gins at the preshock outer boundary Xi where the gas is 
assumed to be unperturbed and the number densities of 
free electrons rio and bound atomic levels rii are in sta- 
tistical equilibrium with radiation field of the shock wave. 
Preshock integration is done on the interval [Xi, Xg_i/2], 
where Xg_i/2 is the space coordinate of the cell center 
just ahead the discontinuous jump. Initial conditions for 
the postshock integration are obtained from the Rankine- 
Hugoniot equations relating variables at the cell centers 
Xg_i/2 and Xg_|_i/2. The outer boundary of the postshock 
region Xjv is set behind the recombination zone. When 
integration of the ordinary differential equations on the 
interval [Xg^if2, Xjy] is completed, we solve the transfer 
equation for the whole slab. Below we discuss the steps of 
the global iteration procedure in more detail. 

2.1. Rankine-Hugoniot relations 

The zonal quantities at the cell centers Xg_i/2 and Xg^i/2 
are related by the radiation modified Rankine-Hugoniot 
equations (Marshak 1958) 

pU^eo = m, (3) 
mC/ + Pg + PR = ei, (4) 
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mi^h+-U^j+Fn + UiEn + PR)^e2, (5) 

where E^, and Pr are the radiation energy density, 
radiation flux and radiation pressure, respectively. 



Pfr^ Pa+ Pc = nukTa + UckTe 



(6) 



is the gas pressure, nn and Ta are the number density and 
the temperature of hydrogen atoms both neutral and ion- 
ized, Te is the temperature of free electrons. The specific 
enthalpy is 



h = - — kTa + - — kTc H \ , 

2 p 2 p P P 

where 

L 

£^cx-XH^(l 



(7) 



(8) 



is the excitation energy, L is the number of bound atomic 
levels. El = xh'^o is the ionization energy, xh is the ion- 
ization potential of the hydrogen atom. 

Eqs. (H) - (H) are solved for the compression ratio 



u- 



(9) 



^ - ^ TT4- ' 

p u+ 

where superscripts minus and plus denote for brevity the 
quantities defined at cell centers Xg_i/2 and re- 
spectively. 

Substituting (^) and into (^) and (^) and then com- 
bining two resulting equations we obtain 

-Br] + C = 0, 



(10) 



where 



. ^ (U-f 2Fr -P+ 2E^+Pr: 

5 5 TO op 



(11) 



5=(a,r + (L/-)V5^-^^i±^, (12) 

P op 



The last two terms in (pi]), (12) and (15) take into ac- 
count the contribution of the radiation field which rapidly 
increases with increasing upstream velocity. However, for 
models considered in the present study this contribution 
was found to be rather small. For example, at the up- 
stream velocity oiUi — 70 km s^^ omitting of these terms 
leads to the increase of the compression ratio 77 by nearly 
one percent and to the increase of the postshock temper- 
ature by less than one tenth of the percent. Neverthe- 
less, we retain these terms in Eqs. (^) and (||) in order to 
use the momentum flux Ci and the total energy flux 62 for 
checking the accuracy of the model throughout the slab. 

2.2. Fluid dynamics and rate equations 

The solution vector of the system of ordinary differential 
equations is 

Y ^ {fii,... ,nL,nc,Ea,Ec,U}, (16) 



where Ea — ^n^kTa and E^ = ^n^kTc are the transla- 
tional energies of hydrogen atoms and free electrons per 
unit volume, respectively. Here and below the tilde denotes 
the quantities expressed per unit mass, that is, fii — Ui/ p, 
Ea = Ea/p etc. 

The assumption of the one-dimensional steady flow 
in planar geometry implies that throughout the slab the 
continuity equation is written as pU = to = Co- Thus, the 
system of the fluid dynamics and rate equations is 



dhi 
~dt 



L 



ni PjL 



dfic \ ^ 

-dT = 2^"'^" 

i=l 



• ^ ^ Pu i J 



i=l 



dE, 
~dt 
dE 



-Pa 



dV_ 
~dt 



J, Pq Qc\c Q'm 

at dt 



IdEn 
p dt ' 



(17) 

(18) 

(19) 
(20) 



dU 

(13) ^ 



1 /2 

and At = (Pg/p) is the isothermal sound speed. 

The postshock electron temperature Tj*" is determined 
in the assumption of the adiabatic compression of the elec- 
tron gas: 



(14) 



where 7 = | is the ratio of speciflc heats. The postshock 
temperature of heavy particles is 

n„ '5 n„k \ rj-^ 



--Fr+ ^ 2 E^, -P{, -E+-P+ 
5 n^kU^ 5 



(15) 



1^ 

to dt 



1 

pc 



XyF^ dv, 



(21) 



where Py and Pji are the total (coUisional plus radiative) 
rates of bound-bound upward and downward transitions 
between the i-th and j-th atomic levels, respectively, P^^ 
and Pfti are the total rates of ionizations and recombina- 
tions for the i~th atomic level, Qe\c and Qinc are the rates 
of energy gain by electron gas in elastic and inelastic col- 
lisions, respectively, is the monochromatic extinction 
coefficient, V — 1/p is the specific volume, c is the velocity 
of hght. 

Eqs. ( p7| ) - (|2^) are written for time derivatives in their 
left -hand-sides, whereas for calculation of the shock wave 
spatial structure it is more convenient to integrate the dif- 
ferential equations with respect to the space coordinate X. 
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Furthermore, the time derivatives of the gas pressure and 
specific volume should be expressed in terms of integrated 
variables. Differentiating ra) and (0) we have 



dPg 2 / dEa dEc 

~dr ~ ~dr ^ IT 



dV_ 
It' 



and 

dV _ J_ /dPg dPn\ 

~dt ydT ^ ~dr) 

Substituting ( p^ ) and 

oc 

dPn U 



(23) 



dt 



X^F^ dv. 



(24) 







into (23) we write the time derivative of the specific vol- 
ume as 



dV _ ,^1 dEg dE, 
dt \ dt dt 



where 



A 



-mU (1 - 



m (1 — /3^) pc 



Xi^F^diy, 



(25) 



(26) 



(27) 



and /3 = ar/U. 

With dl) and (H) Eqs. (0) - (|l|) can be rewritten 

as 



U 



U 



u 



u 



^dfii 
dX 

dfic 
dX 

dE„. 



njP;jt + -hcPKi -ni^ Pij - n^PtK, (28) 



i=i 



1=1 



I - yiPg 

dE, ^ l-APg 
dX I - AP„ 



Q'n\c ^" Qclc L''^9^e 



dU 
dX 

where the terms 
1-AP 



^ dEa „ dEc ^ 
Am—-Am—+<R^, 



(29) 

(30) 
(31) 
(32) 



APa \ J 

I - yiPg 1 - yiPg y u 

APa 1 dEn 



l-AP„rn dt ' 



(33) 



9^0 



APq 



1-AP. 



-Pa 



I -APa 

I - AP^ 



-Pr 



u 



1 - ylPa 1 dEB 



1 -yip. 



dt 



(22) and 



(34) 



(35) 



take into account the contribution of the radiation field 
and can be omitted if i^p Ea + E^ and Pr ^ Pg. 

The excitation rate per atom in the i-th initial state 

is 



Pij — ^cCij 



Ri 



(36) 



and the de-excitation rate per atom in the j'-th initial 
states is 



P — 



(37) 



where n* and n* are the bound-level number densities 
of hydrogen atoms in z-th and j-th state given by the 
Saha-Boltzmann equation for the actual non-equilibrium 
number density of free electrons ric. 
The collisional excitation rate is 



J2 ^iiji' 



exp 



XH 
,■2 



XH 



(38) 



where gj is the statistical weight of level j and Tujii is 
the effective collision strength. In the present study the 
effective collision strengths were evaluated from analytic 
fits by Scholz & Walters (1991) for Is — 2 transitions, 
Callaway & Unnikrishnan (1993) for Is — 3 transitions and 
Callaway (1994) for 2 — 3 transitions. For other transitions 
the effective collision strengths were calculated using the 
Chebyshev polynomial fits of the data by Aggarwal et al. 
(1991). 

The photoexcitation rates are given by 



Rij = in 



hu 



J^dv, 



(39) 



where aijiy) is the spectral line absorption cross-section 
in transition i j and Ji, is the monochromatic mean in- 
tensity. The rate of downward (spontaneous plus induced) 
radiative transitions j ^ i is 



-Rl 



-i47r 



hv 



Ju exp 



hv 



(40) 



The total rate of ionizations from the i-th level is 
P,^ = n^C,^ + (41) 
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and the total rate of recombinations to the i-th level is 
P.i = §^(n,Q, + Rl), (42) 

where n* is the number density of hydrogen ions eval- 
uated from the Saha-Boltzmann equation for the actual 
non-equilibrium number density of free electrons rie- The 
coUisional ionization rate is 



(43) 



where Xi is the potential of ionization from the i-th level, 
the r(Te) is a slowly-varying function of Tc evaluated 
with fitting formulae given by Mihalas (1967). The rates 
of phoionizations and photorecombinations are 



4. / '^^J.d. 
hu 



(44) 



and 



-Rl 



oo 

<4./ 

K J 



ttbf (*) ( '^hv 



hv 



(45) 



where abf(*, J^) is the bound-free absorption cross section 
at frequency by hydrogen atoms in the i-th state. As 
is shown below, the role of coUisional processes in the 
both bound-bound and bound-free transitions through- 
out the radiative shock wave is negligible since their rates 
are smaller by several orders of magnitude in comparison 
with rates of radiative transitions. 

The rate of energy gain by electrons in inelastic colli- 
sions is (Murty 1971) 



^inc — 



dEi dE^^ 1 



dt 



dt 



--V-Fr, 
P 



where 



V- Fr = 47r / (771, - K^J^) dv 



(46) 



(47) 



is the divergence of radiative flux, rju and the emis- 

sion and absorption coefficients, respectively. 

Behind the discontinuous jump the temperature of 
heavy particles exceeds the electron temperature Te 
and free electrons acqiiirc the energy in elastic collisions 
with heavy particles. The rate of energy gain in elastic 
collisions is a sum of the rates corresponding to collisions 
with hydrogen ions and neutral hydrogen atoms: 



— ^ei ~t~ Qe 



(48) 



The rate of energy gain in elastic collisions with hydrogen 
ions is (Spitzer & Harm 1953) 



3/9 4 



where 



252Te 



3/2 



(49) 



(50) 



nelnA 

is the equipartition time and 

A = 9.43 1.15 log (T^/ric 

The rate of energy gain by electrons in elastic collisions 
with neutral hydrogen atoms is 

riH mu 



Qea = ni- 



(51) 



(52) 



where the elastic scattering cross section is (Narita 1973) 

{(JeaV^) = j CTeaV^f{v) dv = 



me 



3/2 



4 + 



24 



{1 + 2- 



(53) 



2.3. Radiation transfer equation 

The quasi static radiation transfer equation for the planar 

geometry is 

dI^l„ 



■^fj,!/ Si/ , 



(54) 



dTv 

where I^,,^ is the specific intensity of radiation at the direc- 
tional cosine jjL = cos(9) at frequency v, S^, = rj^/n^ is the 
monochromatic source function and is the monochro- 
matic optical depth defined as rfr^ = Xi^dX. The abso- 
prtion and extinction coefficient are 



= «bf(i') + K;ff(j^) + K^iiv) 



and 



(55) 



(56) 



where /tbf(z^) is the bound-free absorption coefficient, 
ns{y) is the free-free absorption coefficient, o"R(t') is the 
cross section of the Rayleigh scattering by hydrogen atoms 
in the ground state, ctt is the Thomson scattering cross 
section, Ke,(v) is the spectral line absorption coefiicient. 
Expressions for Ky^iiv) and K^iv) can be found, for exam- 
ple, in the book by Mihalas (1978). The Rayleigh scat- 
tering cross section UYdy) was evaluated using the fitting 
formulae given by Kurucz (1970). 

The continuum emission coefficient is 



2hu^ 



■ L 



(57) 
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where 



4e6 / 27: V'^ gs{v,T,) 



(58) 



is the free-free absorption cross section and gQ{v,Tc) is 
the free-free Gaunt factor. 

The spectral Une absorption and emission coefficients 
are computed in the assumption of the complete redistri- 
bution, that is, 



47r 



1 



nigj 



and 



(59) 



(60) 



where gi and gj are statistical weights, Aji and Bij are 
Einstein coefficients, (jji, is the spectral line absorption 
profile. Thus, we take into account the frequency depen- 
dence of the source function ^'^(i^) ~ rji{v) / ki{u) within 
the spectral line profile. 

The radiative transfer problem is subject to the bound- 
ary conditions implying that the incident radiation on 
both faces of the slab is due to the thermal equilibrium 
radiation from the ambient unperturbed medium: 



I 



fj.uN+1/2 



(61) 



where B^[Ti) is the Planck hmction, Ti is the temperature 
of the unperturbed gas and superscripts minus and plus 
correspond to /i < and /i > 0, respectively. The transfer 
equation ( p4[ ) with boundary conditions (^) is solved for 
the symmetric intensity average 



u 



1/2) 



(62) 



defined in the angle range < /i < 1 (Feautrier 1964). The 
main difficulty arising in the solution of the transfer equa- 
tion for the shock wave is due to the very small optical 
depth intervals Ar^ in the vicinity of = at frequen- 
cies lower than the ionization threshold of the level i — i. 
In the present study this difficulty was overcome with use 
of the improved Feautrier solution proposed by Rybicki & 
Hummer (1991) which provides the better numerical con- 
ditioning and prevents the lost of machine accuracy in the 
matrix elimination scheme. 

In each cell center the intensity average Wp,yj-i/2 is 
determined for sets of frequency and angle points. In or- 
der to compute the frequency integrals the continuum 
frequency range is divided into intervals with bound- 
aries at the threshold ionization frequencies Ui. The upper 
boundary of the continuum frequency range is z/g . Contin- 
uum intervals are represented by nodes of the Gauss- 
Legendre quadratures and spectral lines are represented 
by Tii nodes. The radiation energy density E^, radiative 



flux and radiation pressure Pr are determined by sum- 
mation of frequency integrals over continuum intervals and 
over spectral lines: 



Er^Y.^ I J^dv + \^^ J J^dv, (63) 



Pr = ^ 47r / H^dv + ^ 47r / H^dv, (64) 
— 1 ^ ^ — 1 ^ 



Pr = Y.^ J K,dy + Y.^ J K,dv, (65) 



Ui—/\Ui 



where is the number of frequency intervals, is the 
number of spectral lines, Ai/^ is the half-width of the i- 
th spectral line frequency interval. The moments of the 
radiation field in Eddington's notation are 



Ju = I u^„dfi, 



1 



where 



(66) 



(67) 



(68) 



(69) 



is the antisymmetric flux-like average. 
3. Properties of the solution 

The thickness of the adiabatic compression region where 
thermodynamic variables undergo an abrupt change is 
only a few particle mean free paths. Therefore, in com- 
parison with the whole shock wave the adiabatic com- 
pression region can be treated to a high degree of approx- 
imation as an infinitesimaly thin discontinuous jump. The 
assumption of the steady-state flow allows us to represent 
the shock wave by the flat flnite slab comoving with the 
discontinuous jump, so that the integration of the fluid 
dynamics and rate equations follows the element of gas 
which passes through the slab. The axis of the space coor- 
dinate X is chosen in such a way that the velocity of the 
gas element isU = dX/dt > and the radiative flux emer- 
gent from the preshock and postshock outer boundaries is 
Fi < and F^r > 0, respectively. 
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For all models the space coordinate of the preshock 
outer boundary is Xi — —10^ cm, that is the size of the 
preshock region exceeds by an order of magnitude the 
size of the radiative precursor. The gas element crosses 
the preshock outer boundary at temperature Ti, density 
pi and velocity Ui. The radiative field at the preshock 
outer boundary is non-equilibrium, so that the number 
densities of hydrogen atoms in z-th state n, and the 
number density of free electrons should be evaluated 
from the solution of the equations of statistical equilib- 
rium rather than in LTE approximation. To this end we 
solve Eqs. (|7|) and ([l8| ) with time derivatives equated 
to zero using the Newton-Raphson iteration procedure 
(see, e.g., Mihalas 1978). It should be noted, however, 
that the number densities obtained from the solution 
of the equations of statistical equilibrium do not differ 
substantially from their counterparts given by the Saha- 
Boltzmann equation. For example, the hydrogen ioniza- 
tion degree at the preshock outer boundary is in the range 
8.9-10"'^ < xh < 2.5- 10"^ for upstream gas flow velocities 
15 km s^^ < Ui < 70 km s^^, whereas the LTE hydrogen 
ionization degree is xh = 8.4- 10"'^. Departure coefficients 



(70) 



of the first three levels change from b — {0.9, 1.1, 1.1} to 
b = {0.1,0.3,0.9} for the same range of Ui. 

Of great importance is the choice of the space coor- 
dinate Xjv of the postshock outer boundary. In particu- 
lar, the slab should involve the layers of the full hydro- 
gen recombination since otherwize the total radiative flux 
emergent from the shock is underestimated. Test calcula- 
tions show that for the postshock outer boundary located 
behind the maximum of the hydrogen ionization the to- 
tal radiative flux asymptotically approaches its limiting 
value with increasing Xf^. At the same time, the increase 
of Xjv is accompanied by the growth of the oscillation am- 
plitude of iterated variables in the outermost layers of the 
postshock region, so that for too large X^ the global iter- 
ations ultimately diverge. This occurs mostly due to the 
limited accuracy of integration of the ordinary differential 
Eqs. ( p8| ) - (p^). In the present study the space coordi- 
nate Xj\r was determined for each shock wave model from 
trial computations and ranged from X^ = 5 • 10^ cm for 
Ui = 70 km to Xn = 6- lO'^ cm for Ui ^ 15 km s^^. 

In order to solve the radiation transfer equation ( ^4[ ) 
with the Feautrier technique the slab is divided into 
4000 < < 6000 cells with 1500 cells in the preshock 
region. The cells are smallest at Xg — and increase in 
both directions from the discontinuous jump according to 
the geometrical progression. For all models the cell size at 
the discontinuous jump is AXg_i/2 = ^^3+1/2 = 0-1 cm. 
The continuum frequency range is divided into Lc = L 
frequency intervals [vi, where Vi is the threshold fre- 

quency for ionization from the i~th state, 1 < i < L and 
vq = 10^^ Hz is the the upper boundary of the frequency 



range. The frequency integrals (|6^) - (|65|) were computed 
using the Gauss-Legendre quadratures with ric = 12 
nodes for the continuum intervals and with np = \l nodes 
for spectral lines. For the angular integration of the Ed- 
dington moments (|6l) - (|68|) we make use of the Gauss- 



Legendre quadrature with = 12 nodes. 

In order to consider the dependence of the solution on 
the number of bound atomic levels L and on the width of 
the frequency range [ul, uq] we calculated three sequences 
of the models with L = Lc = 2, 3 and 4. The unperturbed 
temperature and density of the ambient hydrogen gas are 
Ti = 6000K and pi — 10^^" gmcm^"^, respectively. These 
conditions are typical for atmospheres of Cepheids and 
RR Lyr type variables. The convergence of global itera- 
tions is very sensitive to the initial approximation, so that 
the flrst model of each sequence with Ui = 15 km s~^ (the 
upstream adiabatic Mach number is Mi « 1.6) was com- 
puted using the initial LTE approximation with postshock 
temperatures Ta and exponentially decaying with in- 
creasing distance from the discontinuous jump over the 
e-folding distance of 10^ cm. Each converged model was 
used as an initial approximation for the following model 
of higher amplitude. In order to avoid strong initial os- 
cillations of iterated quantites the models were computed 
with amplitude increment AUi — 1 km s^^. The conver- 
gence of global iterations depends also on the size of cells 
approximating the postshock region, so that is why the 
shock models were represented by several thousand cells. 

Hydrodynamic variables that are most sensitive to the 
accuracy of the solution are the electron temperature To 
and the temperature of heavy particles Ta in the vicinity of 
the postshock outer boundary as well as the radiation flux 
-Fr emergent from the slab. Therefore, the convergence of 
the solution on the (fc -I- l)-th iteration was controlled 
by the maximum correction of the electron temperature 
within the slab 



yi = max 



1-T 



(fe) 



(fe+i) 



cj-l/2/-'cj-l/2 



(71) 



and by the correction of the radiative flux Fi emergent 
from the outer boundary of the preshock region 



y2 



(fe) ,p(fe+i) 



(72) 



If global iterations converge, both yi and 2/2 exhibit the 
exponential decay with increasing number of global iter- 
ations k. The typical behaviour of yi and y2 is displayed 
in Fig. 1^. After a few hundred iterations these quantities 
cease to decay and oscillate around the constant average 
values due to the limited accuracy of computations. In the 
present study the stability of the converged solution was 
checked for some models within several thousand itera- 
tions. 

For controlling the degree of approximation we evalu- 
ated also in all cell centers of the model the momentum 
flux Ci and the total energy flux C2 defined by Eqs. ^ 
and (|^), respectively. By definition, both these quantities 
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200 800 
iteration 



1000 



Fig. 1. Convergence plots as a function of the number of 
global iterations for the model with Ui — 59 km s~^. Up- 
per panel: the maximum corrections of the radiative flux 
emergent from the slab (solid line) and the electron tem- 
perature (dashed line). Lower panel: the maximum rela- 
tive change of Ci (solid line) and 62 (dashed line) through- 
out the slab. 

should remain constant throughout the slab if the solu- 
tion is exact. On the lower panel of Fig. |l| are shown the 
maximum deviations 



max 

l<j<N+l 



(« = 1,2), 



(73) 



where Gn are defined at the preshock outer boundary j — 
1 and Cy „i/2 are defined at the cell centers 1 < j < N+1, 
respectively. As is seen from the lower panel of Fig. |l|, the 
total energy flux C2 is constant throughout the slab with 
errors of less than a few tenths of a percent. 



4. General description of models 

Shock wave models computed with the same upstream gas 
temperature Ti and gas density pi can be instructively 
represented by the Hugoniot curve (see, e.g., Zeldovich & 
Raizer 1966) relating the total pressure P and the specific 
volume V ~ 1/p just ahead and just behind the discon- 
tinuous jump. In our case such a representation can be 
only approximate because both the gas density and the 
total pressure ahead the discontinuous jump grow with 
increasing upstream velocity due to radiative heating in 
the preshock region. In particular, for upstream velocities 
15 km s^^ < Ui < 70 km s^^ the gas density increases 
by as much as two percent, whereas the total pressure 
increases by nearly a factor of three. It should be noted 




Fig. 2. The pressure ratio / P versus the inverse 
compression = p~ j for models with L = 3 and 
15 km s-i < [/i < 70 km s^^. 



that the growth of the preshock gas density and the gas 
pressure becomes substantial only at upstream velocities 
U\ > 50 km s~^. In comparison with the amplidute of 
the jump these changes, however, are enough small. The 
Hugoniot diagram for the sequence of models with L = 3 
is shown in Fig. ^. The state of the gas at the cell center 
Xg_i/2 is given for all models by (1,1) and is not shown 
on the plot. The state of the gas just behind the discon- 
tinuous jump at the cell center is represented for 
each model by the filled circle. 

General properties of converged models computed for 
i = 3 and upstream velocities Ui = 20, 40, 60 km s~^ 
are displayed in Figs. where the electron tempera- 

ture To, the temperature of heavy particles T^, the ratio 
of the number density of hydrogen atoms in the second 
state to the number density of hydrogen atoms n2/n}i, 
the hydrogen ionization degree xh, the compression ratio 
p/ Pi, and the total radiative flux Tr are shown as a func- 
tion of the distance from the discontinuous jump X. For 
the sake of convenience the independent variable X is in 
logarithmic scale and each plotted variable is represented 
by two branches corresponding to the preshock and post- 
shock regions, respectively. Just behind the discontinuous 
jump the temperature of neutral hydrogen atoms and hy- 
drogen ions Ta exceeds the temperature of the electron 
gas To. That is why in the upper panel of Figs. |3|-|| the 
postshock temperature plots are represented by two con- 
verging curves. 

Properties of shock wave models are listed also in Ta- 
ble |l|. The column Ui give the upstream velocity at the 
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Fig. 3. The electron temperature Tg, the temperature of 
heavy particles Ta, the ratio of the population number 
density of the second bound level to the number density of 
hydrogen atoms n2 /nn , the degree of hydrogen ionization 
xa, the compression ratio p/pi, and the total radiative 
flux Fft versus the distance from the discontinuous jump 
X for the model with L = 3, Ui = 20 km s~^. In dashed 
and dot-dashed lines are given the plots for the preshock 
and postshock regions, respectively. 



outer boundary of the slab, Fji is the radiative flux emer- 
gent from the preshock outer boundary and Ftx/G2 is the 
ratio of the radiative flux to the total energy flux at the 



preshock outer boundary. Columns headed by X 



a-i/2 



de- 



scribe the state of gas just ahead the discontinuous jump. 
Because in the preshock region the divergence of radiative 
flux is negative, we use its absolute value T> — |V- F^]/p. 
The total rise of the gas temperature and the total rise 
of the hydrogen ionization degree in the preshock region 
within the interval [Xi, Xg_i/2] are given in columns AT 
and Axh, respectively. The columns headed by Xg^i/2 
give the temperature of heavy particles T^, the temper- 



4x10* 



3x10* 



H 2x10' 
lO'' 
10"^ 
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bxlO"* 
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0.4 

K 
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Fig. 4. Same as Fig. || for model with L 
40 km s"^ 



= 3, C/i = 



ature of the electron gas and the compression ratio 
rj = / p~ just behind the discontinuous jump at the cell 
center Xg_^^i/2- The columns headed by maxD give the 
maximum value of the divergence of radiative flux in the 
postshock region (in these layers V- Fr > 0) and the space 
coordinate X where this maximum is attained. The last 
two columns are headed by maxxn and give the maxi- 
mum ionization degree of hydrogen xh in the postshock 
region as well as the distance X of this maximum from 
the discontinuous jump. 

A cursory inspection of Table |l] shows that the total 
radiative flux Fr nearly does not depend on the number of 
bound atomic levels L. In other words, the radiative flux 
emergent from the shock wave is independent of the lower 
limit VL of the frequency range provided that vl < V2- 
At the same time, the change of vl is accompanied by re- 
distribution of the energy of radiation field within [v^ , t'o] 
which is most perceptible within the Lyman continuum 
In particular, for the narrower frequency range (e.g. for 
L = 2) we have both stronger radiative heating and higher 
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Table 1. Parameters of shock wave models with pi = 10 gmcm ^ and Ti = 6000K. 
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hydrogen ionization in the preshock region. The structure 
of the postshock region is also sensitive to the width of 
the frequency range [j/LjJ/o]. In particular, for the nar- 
rower frequency range the thermal equipartition and the 
maximum of ionization are attained at larger distances 
from the discontinuous jump. 



5. Radiation field 

We consider the shock waves propagating through the par- 
tially ionized hydrogen gas which is opaque both in the 
Lyman continuum and in the Lyman lines. The total op- 
tical depth of the slab at the frequency v\ corresponding 
to the threshold ionization from the ground state is in the 
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Fig. 5. Same as Fig. || for model with L = 3, Ui = 
60 km s"^ 



range 10^ < '''{'^i) ^ 10^, the maximum t{i'i) ^ 10^ being 
reached at Ui f« 40 km s~^. At the same time the slab is 
nearly transparent for continuum radiation at frequencies 
v < vi. The optical depths in the centers of Balmer lines 
are ^ 1 in the preshock region and gradually increase 
with increasing distance from the discontinuous jump in 
the postshock region. At distances X > cm the post- 
shock optical depth in the center of the Ha spectral line 
measured from the discontinuous jump is T,y > 1. In Fig. || 
we give the plots of monochromatic optical depths mea- 
sured from the discontinuous jump to the layer with space 
coordinate X for the heads of continua with threshold fre- 
quencies i^i (1 < j < 3) and for the centers of spectral lines 
Lya and Ha. 

Thus, the Eddington approximation cannot be applied 
for the solution of the radiation transfer problem in shock 
waves. Indeed, it is valid only within the Lyman contin- 
uum whereas the most of radiation is transported at fre- 
quencies lower than the Lyman edge where the angular 
dependence of the specific intensity I^^ must be taken 
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Fig. 6. Monocromatic optical depths t^, measured from 
the discontinuous jump for the preshock (upper panel) and 
postshock (lower panel) regions of the model with L = 3, 
U\ — 40 km s^^. In solid lines are shown the optical depths 
at the heads of the continua and in dashed lines are shown 
the optical depths at the centers of spectral lines. 

into account. Moreover, the correct calculation of the an- 
gular dependence of the specific intensity I^^y is of great 
importance because for v <v\ the term 

- f S,it)e-^'-^"^^f'dt (74) 



rapidly decreases for /i ^ |1| and the Eddington factor is 
U = K,/J, < 1/3. 

The interaction between the radiation field and the 
gas material is described in terms of the divergence of 
radiative flux V-Fr. By definition, this quantity is negative 
when the gas absorbs more energy than it emits, whereas 
when V • Fr > 0, the gas radiatively cools. The typical 
plot of the divergence of radiative flux expressed per unit 
mass i V- Fr is shown on the upper panel of Fig. ^ for the 

shock wave model with L ~ 3 and Ui = 40 km s^^. The 
gradual decrease of ^ V-Fr in the preshock region is due to 
the strong absorption of the Lyman continuum radiation 
which leads to the growth of the gas temperature when 
the parcel of gas approaches the discontinuous jump. The 
heating of gas rapidly increases with increasing upstream 
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Fig. 7. The divergence of radiative flux per unit mass (up- 
per panel) and the radiation energy density (lower panel) 
as a function of the distance from the discontinuous jump 
for the model with L = 3, J7i = 40 km s~^. In dashed and 
dot-dashed lines are given the preshock and postshock 
plots, respectively. 



velocity and as is seen from Table 0the quantity 2) = |V- 
FrI/p at Xg_i/2 changes be several orders of magnitude 
for 15 km s'^ < C/i < 70 km s"^ 

The jump of ^V- Fr within [Xg_i/2, -^^3+1/2] is due 
to the adiabatic heating of electrons at the discontinu- 
ous jump. The layers of the most efficient radiative cool- 
ing in the postshock region with fastest decrease of the 
gas temperature are revealed by the maximum of ^V - Fr 
whicj also very rapidly increases with increasing Ui (see 
Table |). The diver gence of radiative flux approaches zero 
with increasing distance from the discontinuous jump in 
the both preshock and postshock regions. This implies 
that at large distances from the discontinuous jump the 
total radiative flux tends to be constant. 

While the parcel of gas approaches the discontinuous 
jump in the preshock region the radiation energy density 

gradually increases due to the radiative heating of gas. 
Behind the discontinuous jump E^i continues to increase 
and reaches the maximum in the layers of the fastest re- 
combination of hydrogen atoms. The maximum of the ra- 
diation energy density Sr. in the postshock region approx- 
imately corresponds to the point where the total radiative 
flux is Fr = 0. It should be noted, however, that this con- 
dition is exact only in the diffusion approximation (see, 
e.g., Mihalas & Mihalas 1984) and as is shown below the 
space coordinate of the layer with = depends on the 




40 50 
Uj (km s"') 

Fig. 8. The ratio of the radiative flux Fr emergent from 
the preshock outer boundary to the total energy of the 
shock wave C2 versus the upstream gas flow velocity Ui 
for models with L — 3. Along the upper horizontal axis 
are given the upstream adiabatic Mach numbers Mi. 



frequency ly. The typical plot of Fr is shown on the lower 
panel of Fig. 0. For all models considered in the present 
study the radiation energy density Fr can be neglected 
in comparison with the thermal energy of the gas. For ex- 
ample, in the model represented in Fig. Q the ratio of the 
radiation energy density to the enthalpty is Fr//i < 10~^ 
throughout the whole shock wave. 

The radiative flux outgoing from the preshock outer 
boundary is nearly the same as that outgoing from the 
postshock outer boundary, that is, — Fi « F/y. The radia- 
tive flux emergent from the shock wave rapidly increases 
with increasing upstream velocity and at Ui ~ 60 km 
more than 90 percent of the total energy flux is contained 
in the radiation flux. As is seen from Fig. |^ where we give 
the plot of the ratio Fr/C2 versus the upstream velocity 
Ui, at upstream Mach numbers AIi > 6 the major part of 
the shock wave energy (i.e. more than 90%) is lost due to 
radiation. 

At upstream velocities Ui < 40 km s^^ the total ra- 
diative flux is nearly constant in the preshock region and 
undergoes the substantial changes only behind the discon- 
tinuous jump in the layers where the gas radiatively cools 
and recombines. However the role of the Lyman contin- 
uum radiation rapidly increases with increasing upstream 
velocity and for C/i > 40 km s^^ the total radiative flux 
Fr shows perceptible changes in the preshock region (see 
Fig. I). 

Of great interest are the spectral distributions of 
monocromatic quantities describing the radiation field as 
a function of frequency v. In Fig. ^ we give the spectra 
of monocromatic radiative flux F^ for three different lay- 
ers: the preshock outer boundary Xi, the discontinuous 
jump Xg and the postshock outer boundary Xn in the 
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Fig. 9. Monochromatic radiative flux F^, as a function of 
frequency v for the preshock outer boundary Xi, discon- 
tinuous jump Xg and postshock outer boundary . For 
the sake of convenience the plots are arbitrarily shifted 
along the vertical axis. 

model with L — Z and U\ = 40 km s^^. As is seen from 
these plots, the most of radiation is transported within 
the frequency interval vi < v < v\ limited by the edges 
of the Balmer and Lyman continua. The contribution of 
the Lyman continuum radiation is perceptible only in the 
vicinity of the discontinuous jump. 

The conspicuous feature of these plots is also that the 
Ha spectral line is in emission throughout the shock wave, 
whereas the Lya and Ly/3 lines are mostly in absoprtion 
and are revealed in emission only in the vicinity of the dis- 
continuous jump. Within the range of upstream velocities 
15 km s~^ < C/i < 70 km s~^ the frequency dependen- 
cies of the radiative flux are qualitatively similar to that 
shown in Fig. |^. 

The role of various spectral intervals in the total radi- 
ation transfer can be evaluated with use of terms 

Vi-X 

F^^ = 'in J H^dv (75) 
and 

= At: j H^dv (76) 

representing the radiative flux integrated within the con- 
tinuum frequency intervals [vt^vi^i] and within the spec- 
tral lines, respectively. Here /S.Vi is the half-width of the 



L=3 Uj = 20 km s' 




- 



2 4 6 

log |X| 

Fig. 10. The radiative flux in continuum intervals F^ (up- 
per panel) and in spectral lines (lower panel) for model 
with L = 3, ?7i = 20 km s~^. In dashed and dot-dashed 
lines are represented the preshock and postshock plots, 
respectively. 

spectral line frequency interval. Figs. |l^-|l^ show Ff and 
-F"/ as a function of the distance from the discontinuous 
jump for models with Ui = 20, 40 and 60 km s~^, respec- 
tively. For the sake of convenience of the graphical rep- 
resentation the fluxes are given in the logarithmic scale. 
As is seen from these plots throughout the shock wave 
the most of radiation is transported within the Balmer 
continuum. The role of the Lyman continuum becomes 
perceptible at upstream velocities Ui > 40 km s^^ only 
in the vicinity of the discontinuous jump: — 10** cm < 
X < 10'* cm. In the preshock region the radiative flux 
is constant at frequencies v < i^i. The Lyman continuum 
flux {vi <v< vq) rapidly increases with approaching the 
discontinuous jump and becomes perceptible at distances 
smaller than X ^ 10** cm. But even in the vicinity of 
the discontinuous jump the radiative energy transport in 
the Lyman continuum remains small in comparison with 
that at lower frequencies until the upstream velocity is 
t/i < 40 km s-i. 

As is seen from the plots given on the lower panels of 
Figs. |o|-|l^ the contribution of the spectral line radiation 
into the total radiative flux is negligible. Throughout the 
slab the most of the spectral line radiation is transported 
by Ha and only at upstream velocities /7i > 50 km s^^ in 
the close vicinity of the discontinuous jump the spectral 
line radiation is transported mostly in Lya. 
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Fig. 11. Same as Fig. |10| for model with L = 3, Ui = Fig. 12. Same as Fig. |10| for model with L = 3, Ui = 
40kms-i. 60kms-i. 



Behind the discontinuous jump the monochromatic ra- 
diative flux undergoes substantial changes within the ra- 
diative cooling region from where the most of radiation 
is emerged in both opposite directions. By definition, the 
monochromatic radiative flux changes its sign in the layer 
where V • Fi, reaches the maximum. In plots given in 
Figs. 10 -|l2| these layers are revealed as deep maximuma 
of the continuum and spectral line fluxes F[ and Ff, re- 
spectively. 

As is seen, the region emerging radiation is spread 
along the space coordinate in wide ranges: 3.8 < logX < 
6.5 for Ui = 20 km s'^ and 3.5 < logX < 4.1 for 
J7i = 60 km s^^. This is obviously due to strong frequency 
dependence of optical depths in layers emitting the ra- 
diation. 

6. Preshock region 

Effects of the electron thermal conduction in the vicin- 
ity of the discontinuous jump are not considered in the 
present study, therefore throughout the preshock region 
the translational energies of the both hydrogen atoms and 
free electrons are equal and can be described in terms of 
the gas temperature: T = = T^. For T = 6000K and 
p — 10~^° gm cm^^ most of the hydrogen atoms are in 
the ground state, so that the zone in which the radiation 
field interacts with gas material is only of several mean 
free paths of photons at the frequency of the Lyman con- 
tinuum edge vi. For our models this corresponds to the 
distance from the discontinuous jump of X '--^ 10^ cm. 



Within this zone known as the radiative precursor the 
gas undergoes both radiative heating and photoionization. 
The rize of the gas temperature AT and the hydrogen ion- 
ization degree Axh within the preshock region is nearly 
negligible in shocks with upstream velocities smaller than 
Ui ~ 30 km s~^ but at larger upstream velocities AT 
and Axh very rapidly grow with increasing Ui (see Ta- 
ble |l|). At upstream velocities Ui > 65 km s~^ the hy- 
drogen ionization degree ahead the discontinuous jump is 
xh > 0.5 and the gas material becomes more transparent. 
This accounts for the small decrease of 2) — |V- Fr |/p for 
Ui > 65 km s-^ 

Throughout the preshock region the coUisional ioniza- 
tions can be neglected since their rates are by several or- 
ders of magnitude smaller than those of photoionizations. 
In Fig. |l^ we give the plots of the net rates of radiative 
ionizations 



(77) 



for first three atomic levels in models with Ui = 20, 40 
and 60 km s^^. At distances from the discontinuous jump 
less than X ~ 10^ cm the absorption coefficient in Lya 
and Ly/3 lines exceeds that of the Lyman continuum. For 
upstream velocities Ui < 30 km the ionization from 
the ground state is still small and in layers with X < 
10^ cm the hydrogen ionizes mostly from excited levels 
i > 2. However the rapid growth of the ionization from 
the ground state with increasing upstream velocity cancels 
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Fig. 13. Net radiative ionization rates for the ground 
state (solid line), second level (dashed line) and third 
level (dot-dashed line) in the preshock regions of the 
shock wave models with Ui = 20 km s^^ (upper panel), 
Ui = 40 km s'^ (middle panel) and Ui — 60 km s^"'^ 
(lower panel). 



this effect and at t/i > 30 km s ^ the hydrogen ionizes 
mostly from the ground state. 

The columns AT and Axh of Tabic |l] give the total 
growth of the gas temperature and ionization degree in the 
preshock region. As is seen, at constant upstream velocity 
Ui both AT and Axh decrease for larger L. As was noted 
above the total radiative flux emitted by the shock wave 
nearly does not depend on the frequency range [vl,vq\ 
for L > 2 and for smaller the monocromatic flux is 
redistributed among lower frequencies. In particular, the 
fraction of the Lyman continuum radiation in the total 
radiative flux decreases. 

The heating of the gas within the radiative precursor 
is accompanied by the rise of the gas pressure Pg and the 
gas density p. However the increase of p is not large and 
for Ui = 60 km s~^ the gas density just ahead the dis- 
continuous jump is p~ = 1.0085pi, so that the upstream 
gas flow velocity decreases by less than one per cent and 
is U~ = 0.9916J7i. More important is the rise of the gas 
temperature T because it leads to the slow decrease of the 
compression ratio p'^ / p^ with increasing Ui in shocks with 
upstream velocities Ui > 55 km (see Table |^). Such 



a decrease of the compression ratio p'^ / pi with increas- 
ing upstream velocity is easily explained in terms of the 
Rankine-Hugoniot relations (^ - (||)- Adopting Fj^ = 



in (11) and neglecting the contribution of the radiation 



energy density and radiation pressure Pr in Eqs. (11) 
and (|l2|) we find that the compression ratio at the discon- 
tinuous jump is 



4([/ 



-\2 



5{a^y + {U-Y 



(78) 



Thus, the slow decrease of the compression ratio p^/p^ 
with increasing Ui is due to the rize of the isothermal 
sound speed ahead the discontinuity jump. Inspection of 
Table 1^ shows that the effect of decreasing becomes 
less prominent in the models computed with larger L and, 
therefore, having somewhat weaker radiative heating of 
the gas just ahead the discontinuous jump. 

For gas densities of 10^ gm cm^'^ S P S 
10^^ gm cm~'^ the dominant recombination mechanism 
is the three-body recombination (see, e.g., Zeldovich & 
Raizer 1967) and as was shown by Gillet & Lafon (1990) 
the precursor heating is due to the partial transformation 
of the ionization energy into the thermal energy via the 
three-body recombination processes. Kuznetsov & Raiser 
(1965) have shown that electrons are captured on the level 
with energy Ui = xh/*^ close to the average kinetic en- 
ergy of free electrons kTe. Thus, i = ^/xHjkTc and for 
the precursor temperature range 6000K < < 10'*K (see 
Table |l|) the free electrons should be captured onto levels 
4 < i < 5. Thus, the atomic model represented by 4 bound 
levels seems do not underestimate significantly the rates of 
three-body recombinations and, therefore, the precursor 
heating is treated quite correctly. This conclusion is con- 
sistent with decrease of the temperature rize AT obtained 
for atomic models with L = 3 and L = 4 (see Table p. 



7. Postshock region 

Just behind the discontinuous jump the temperature of 
neutral hydrogen atoms and hydrogen ions r+ exceeds 
the temperature of free electrons Tj*" (see Table |l|) and 
as is seen from Figs. ||-|| both these temperatures tend 
to equalize with increasing distance X. Together with ex- 
citation of bound atomic levels the thermal equilibration 
is the fastest relaxation process behind the discontinuous 
jump. For the models considered in the present study the 
hydrogen ionization degree ahead the discontinuous jump 
is > 10~^, so that ionized hydrogen atoms are enough 
numerous behind the discontinous jump and free electrons 
gain the energy from heavy particles mostly in elastic col- 
lisions with hydrogen ions. It should be noted that for all 
models throughout the postshock region Qd S> Qca- 

When the temperatures of heavy particles and free 
electrons equalize, the gas radiatively cools and the gas 



16 



Yu.A. Fadeyev & D. Gillet: The structure of radiative shock waves. II 



temperature T = Ta = gradually decreases with in- 
creasing distance from the discontinuos jump. The radia- 
tive cooling is fastest when the divergence of radiative flux 
( ^ ) reaches its maximum. The maximum postshock val- 
ues of iV- Fr and the space coordinate X, where this 
maximum is reached, are given in Table |. The maxi- 
mum of ^V- Fr is prominent only when the postshock 
ionization degree is xu < 0.8 and for larger cch the spa- 
tial dependence of -!-V- Fr shows the plateau rather than 
the maximum. It s&ould be noted that the maximum of 
V- Fr gives only the frequency averaged location of the 
layers emitting the radiation. As was shown above the 
space coordinate of these layers depends on the frequency 
of radiation v. 

The ionization of hydrogen atoms is the slower relax- 
ation process in comparison with thermal equilibration 
and the maximum of xh is reached behind the maximum 
of iV- Fr (see Table p. Throughout the postshock re- 
gion the rates of coUisional processes are by several orders 
of magnitude smaller than those of radiative processes 
and, therefore, can be neglected. Fig. |lj shows the net 
photoionization rates Zt^, for the first three bound levels 
of the hydrogen atom in models with Ui = 20, 40 and 
60 km s^^. In the postshock region the bound atomic lev- 
els are excited faster than hydrogen ionizes, so that at 
upstream velocities Ui < 50 km s~^ hydrogen atoms ion- 
ize mostly from bound levels i > 2. At upstream velocities 
Ui > 55 km s^^ more than ten percent of hydrogen ionizes 
in the precursor and in the postshock region the hydrogen 
ionizes mostly from the ground state. At C/i > 55 km s~^ 
the maximum postshock ionization degree is xh ~ 1- This 
leads to the slower postshock temperature decrease until 
the hydrogen ionization degree is near the maximum (see 
the upper panel of Fig. |^) . 

The full thickness of the relaxation zone behind the 
discontinuous jump is determined by the slowest relax- 
ation process which in our case is the recombination of 
hydrogen atoms. Substitution of the translational ener- 
gies of hydrogen atoms and free electrons = Ea + E^ 
into the Rankine-Hugoniot equations (^) and (||), where 
we neglect terms with E^ and Pr, gives 



(yEt — E'ti^ + ^i?ex — -Ecxi^ + (^Ei — El 



p m 



(79) 



Here subscript 1 denotes the quantities defined at the 
preshock outer boundary Xi . For strong shock waves with 
Eti < Et, Ec-Ki < i^ox and En < Ei the relation is 
rewritten as 



E, 



Ft 



mEt 



(80) 



The first term in (80) is the limiting compression ratio 



(7 -f l)/(7 — 1) = 4 for the strong adiabatic shock front 
{Ml ^ 1). In fact, because of radiative processes in the 
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Fig. 14. Net radiative ionization rates from the ground 
state (solid line), from the second level (dashed line) and 
from the third level (dot-dashed line) in the postshock 
region of the shock wave models with Ui — 20 km s^^ 
(upper panel), Ui = 40 km s~^ (middle panel) and Ui = 
60 km s^^ (lower panel). 



precursor the compression ratio at the adiabatic shock 
front of our models is as large as j p~ =3.6 (see Ta- 
ble |l|). The second and the third terms in ( |80| ) describe 
the additional gas compression which occurs behind the 
discontinuous jump in the postshock region due to excita- 
tion, ionization and recombination of hydrogen atoms. It 
should be noted that because Jri < 0, the third term in 
( ^0|) is always positive. Thus, the transfer of the part of the 
gas fiow kinetic energy into internal degrees of freedom of 
atoms and into radiation field produces the gas density in- 
crease behind the discontinuous jump which is much larger 
than that at the adiabatic shock front. For example, in the 
model with U\ = 60 km s~^ we obtained the final com- 
pression ratio as high as pj pi ~ 65 (see Fig. ||) which 
approaches the maximum compression ratio of « 84 
in the isothermal shock wave where all the shock energy 
is transformed into the radiation. 
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8. Conclusion 

In the work presented here we have siibstantially improved 
the computer code for calculating the structure of radia- 
tive shock wave models and considered almost the entire 
zone of postshock hydrogen recombination in shock wave 
models with upstream Mach numbers 1.6 < Mi < 7.5. 
This allowed us to obtain the reliable estimates of the ra- 
diative flux emitted by the shock wave as a function of the 
upstream gas flow velocity. The most striking result of our 
calculations is that the ratio of the radiative flux to the 
total energy flux of the shock wave very rapidly enlarges 
with increasing upstream velocity, so that for adiabatic 
Mach numbers Mi > 7 (i.e. at Ui > 65 km s~^) the 
major part of the shock wave energy (more than 90%) is 
irreversibly lost due to dissipation processes. 

Another noteworthy feature of our results is the negli- 
gible role of coUisional processes in the both bound-bound 
and bound-free transitions in comparison with radiative 
transitions throughout the whole shock wave. This circum- 
stance makes our models independent of uncertainties in 
cross sections of collisional processes. 

The results described in the present paper were ob- 
tained for the only case with gas temperature Ti = 6000K 
and gas density pi = 10~^° gmcm""^ of the ambient unper- 
turbed medium which is most typical for atmospheres of 
classical Cepheids and RR Lyr type variables. Obviously, 
many more models have to be constructed for various gas 
temperatures and gas densities in order to more clearly 
understand the properties of radiative shock waves ob- 
served in astrophysical phenomena. Of particular interest 
are the gas temperatures and gas densities typical for pul- 
sating late-type stars possessing strong shock-driven stel- 
lar winds. The efforts of further calculations will be also 
directed to the more accurate modelling of strong shock 
waves in order to consider in detail the run of the ratio 
-Fr/C2 as a function of the upstream velocity at Mach 
numbers Mi > 7. 
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